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O ' Abstract 

Neutrino emissivities in a neutron star are computed for the neutrino bremsstrahlung process. In 
the first part the electro-weak nucleon-nucleon bremsstrahlung is calculated in free space in terms 
of a on-shell T-matrix using a generalized Low energy theorem. In the second part the emissivities 
are calculated in terms of the hadronic polarization at the two-loop level. Various medium effects, 
such as finite particle width, Pauli blocking in the T-matrix are considered. Compared to the 
pioneering work of Friman and Maxwell in terms of (anti-symmetrized) one-pion exchange the 
resulting emissivity is about a factor 4 smaller at saturation density. 
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I. INTRODUCTION 



The cooling of neutron stars proceeds via the weak interaction. Since in general one-body 
processes are kinematically forbidden the dominant reactions are assumed to be the neutral 
current two-particle processes 

n + n — > n + n + Uf + Uf, (1) 



n + p — > n + p + uj + Uf, (2) 

and the charged current "modified URCA" process 

N + n^N + p + u e + e~. (3) 

Standard cooling scenarios are mostly based upon the pioneering work of Friman and 
Maxwell In essence their approach amounts to a convolution of the soft free space 
neutrino pair emission and two-body (modified) URCA processes (PCU) with a finite tem- 
perature free Fermi-gas model using Fermi's golden rule to obtain the emission rate. In 
doing so a number of simplifying assumptions were made; in particular (i) the two-body 
interaction between the nucleons was approximated by a central Landau interaction plus a 
one-pion exchange to represent the tensor force, (ii) and only the non-relativistic limit was 
considered, (iii) since it is based upon the quasi-particle approximation non-perturbative 
effects such as the LPM effect were not taken into account, (iv) other medium effects such 
as Pauli blocking in the strong interaction were neglected. It is the aim of the present paper 
to investigate and possibly improve these assumptions. 

In the first part we consider the reactions (PEJ) in free space. Using the fact that the 
energy release in the bremsstrahlung process is very small we apply the soft bremsstrahlung 
formalism of Hanhart et al. Q] and Timmermans et al. 0. This allows one to express 
the bremsstrahlung process in the soft limit model independently in terms of an on-shell T- 
matrix, i.e. phase shifts. In this way we are able to judge the accuracy of past bremsstrahlung 
calculations, which were mostly based upon the use of a one pion exchange (OPE) approxi- 
mation in the non-relativistic limit In the latter case simplifications occur such as the 
vanishing of the vector current matrix elements. 

In the second part we consider the process in the medium. To describe the cooling 
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process of neutron stars through neutrino emission the application of Fermi's golden rule 
in the quasi particle approximation (QPA) was mostly used in the past. To compute emis- 
sivities beyond QPA one needs to start from quantum transport equations. The essential 
physics is then contained in the neutrino self-energies, which appear in the loss and gain 
terms. We will compare the diagrams at the hadron two loop level. It appears that only in 
lowest order in the imaginary part of the hadronic self-energies the use of closed diagrams 
and the application of Fermi's golden rule coincide. 

From the generalized Low-energy theorem it follows that the use of the QPA leads to 
a infrared divergent amplitude, 1/uj. The latter is predicted to be quenched \^ in a medium 
whenever the mean free path of the nucleons becomes on the order of the formation length 
of the lepton pair. This is also known as the Landau-Pomeranchuk-Migdal(LPM) effect in 
case of electromagnetic interactions). The importance of this effect we study by including a 
finite single particle width (imaginary part of the self-energy) which depends on energy and 
temperature. 

In practice in calculating the collision integral one needs to specify the appropriate dia- 
grams and make assumptions about hadronic interactions. In doing so one must be careful 
that symmetries like gauge invariance of the vector current are not violated. We also esti- 
mate the Pauli blocking by replacing the T-matrix by a in-medium G-matrix. In Sedrakian 
and Dieperink [f| the neutrino emissivity was computed including the LPM effect, however 
in the OPE approximation 

Many properties of superfluid matter such as pairing are still known with large uncer- 
tainty. Therefore only non-superfluid matter will be considered. For recent papers about 
pairing we refer to Gusakov Q and Yakovlev et al. . 

Although we will apply the present formalism to neutrino pair emission in neutral weak 
current processes, it is equally valid for soft electromagnetic bremsstrahlung. 

This paper is organized as follows. In section 2 we discuss electroweak bremsstrahlung in 
free space; in section 3 the in-medium process is discussed at the two-loop level. In section 
4 results are presented showing the effects of various approximations. In the appendix a 
summary of quantum transport theory and finite temperature Green functions is presented. 
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FIG. 1: Diagrams for neutrino pair bremsstrahlung of order 1/uj with radiation from external legs. 
II. ELECTROWEAK BREMSSTRAHLUNG IN FREE SPACE 
A. Soft electroweak bremsstrahlung amplitude 

The vv pair emission in a neutron star is characterized by a very small energy transfer 
(on the order of the temperature T ~ 1 MeV), much smaller than any other scale in the 
process like m, or pp. Therefore it is natural to consider the NN — > NNvV process in the 
ultra-soft limit. For simplicity and also to be consistent with the low density limit of the 
medium we will first consider this process in free space. 

Here the treatment of soft NN electroweak bremsstrahlung, discussed in more detail in 
Re, k is —ized. A„ alog ou, y to the electronic _aMun g (Low, Q), 
the first two terms of the expansion in powers of the energy-momentum transfer \q\ < uo 
of the electroweak bremsstrahlungs amplitude are determined by the amplitude for the 
corresponding non-radiative process M = A/ 'uj + B + 0(u>). In the ultra-soft regime (u/p << 
1, where p is the nucleon momentum) the B and higher order terms can be neglected. The 
amplitude of the diagrams in Fig. ^with radiation from external legs only is given by Q], 

M^' a = T 1 S(p 1 - q)Y a v + YlS{p\ + q)T[ + {1 «- 2}. (4) 



4 



The A term in the Low expansion is obtained by considering the limit \q\ < ui — > of 
u>M^ xt ' a ; to this end we expand the various T's with one nucleon off its mass shell, 

T x =< pi,p' 2 |T|pi - q,p 2 >, T[ =< pi + q,p' 2 \T\p 1 ,p 2 >, (5) 

around the on-shell point To, 

Tx = T - g.^-To + T[ = T + g.^To + ... , (6) 
dpi op\ 

and the nucleon propagator S as 

S(p± ,).^)±^ K± 5^M + 0(l) (7, 
p±q — m 2p.q 

with A±(p) = Zgt. 

The hadronic weak interaction vertex in the limit q — > is given by 

n = ^(cv-c Al5 )^, (8) 

where Gf is the Fermi weak coupling constant and r is the isospin-operator. The vector and 
axial-vector coupling constants cy and ca are for neutrons c v = — 1; c\ = —qa = —1.26 
and for protons Cy = 1 — 4 sin 2 Qw « 0.08; = gr^ = 1.26. 

Since the initial/final particles are on mass-shell one has the relations (jf> + m)j u u(p) = 
2p v u(p) and u{p)'ju{p' + m) = 2p u ll(p), which are useful for the vector current. As a result 
in the ultra-soft region (q/p << 1) the vector and axial- vector current matrix element are 
given by 

M V,a = GFCv(_ ^ T aP^ T \ + {1 ^ 2} (9) 

2V2 V Pl-Q PU 



and 

A 2mG F c A ( A+(pi) A+Cpj) \ , , , , 

M u = 7= - T o^ 7^75^ + 7^75^ -7T-, T +{1^2}, (10) 

respectively. Naturally the vector current is conserved: q v M^> a = 0. 



B. Structure of the elastic NN scattering amplitude 



It is clear that the amplitudes in Eqs. © and 
of T . For the elastic AW scattering amplitude Q 5 



IUJ) depend on the Lorentz structure 
lll | for the process AT(p 1 ) + N{p 2 ) — > 



N(p' 1 ) + N(p' 2 ), the covariant form of the on-shell T-matrix can be expressed as 

T = Tdir + Texch= ^ j2FP(s,t,u)[^ P ' 2 )nMP2mp'i)n a u(pi) 

1=0,1 a=l 

+ (-) a U(p / 1 )fi QM (p 2 )M(p / 2 )fi Q M (j9 1 )] J B / , (11) 

where the five Fermi covariants are 

Q a = (fit, Q 2 , ^3, ^4, ^5) = (1, <W ^ 757^, lv, lb)- (12) 
The projection operators on iso-singlet and iso-triplet states are 

B = (1 - n.^/A, B 1 = (3 + t 1 .t 2 )/A, (13) 

respectively. F^\s, t, u) are the invariant functions of the Mandelstam variables s = — (p± + 
P2) 2 , t = —{p'i — Pi) 2 , and u = —(p' 2 — Pi) 2 - For the nnvv and npvv process the isospin 
combinations needed are 

F<r> (s, t, u) = F^\s, t, u) = F« (s, t, u) 

f {up) = ( F W( s ,t,u) + Fi \s,t,u))/2 (14) 

for a = 1, 5. For later use, it is convenient to put the spinors in the exchange term in the 
"normal order" by introducing the functions 

(s, t, u) = F«(a, + £ {-\yCaf>Fj{8, t, it) , (15) 

(3=1 

where C a @ are elements of the Fierz transformation, the explicit form is given ^0, 11 1. Then 
Eq. (|llj) can be rewritten as 

T= £ E^^^^Wa)^^)^)^^)^- (16) 

7=0,1 a=l 

Since for a comparison we will need the cross section in the non-relativistic limit, we also 
give the required non-relativistic decomposition of T (we will reserve latin indices for the 
non-relativistic T-matrix) 
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T = J2%(s,t,u)O v , (17) 

v=l 

where 

% = {%.■> I2, %, %) = (Tc, Tq, T T i, T T2 , T S o) (18) 



and the five independent two-body operators 

O v = (1, o\ • n o% • ft, o\ • k &2 • k, a\ • k' <r 2 • k', S\ • n + a% • n) (19) 

with k = (pi —Pi)/Wi ~ Pi I) k' = {pi +pi)/|pi + pi| and n = (k' x k)/\k' x fc| in the cm- 
system. The terms T c ,Tq, T T1 and T so corresponds to the central, quadratic spin orbit, 
tensor and spin-orbit force, respectively; further we have a second tensor T?2 (instead of the 
spin-spin force). 



C. The nnvv process 



We first treat the n+n^n+n+v+v process. The vector current amplitude follows 
from Eq. © 



M, 



v 



G F Cy / _ Plu_ P[u_ _ P2v_ P^_ 

2y/2 ^ pi.q p[.q p 2 -q p'-i-q 



£ ta)n a «(p 2 )«(pi)n°tt(pi) + R a {p' 2 <-> p'j' 



(20) 



a=l 



The axial-vector current amplitude follows from Eq. (jl(J|) 

2mG F g A ' 



Aft 



^ E ^ [«(p' 2 )0 QU (p 2 )«(pi) ( - n«^7„ 
2y 2 a= i L ^Pi-9 

L+(pi 

2p[.q 

v 2p v q 2p' 2 .q 



75 



+7,75^^^)«(pi) + (-Wi^fe) 



+ (1 «-> 2) 



A+fa) a+(p 2 ; 



(21) 



For later use we also give the non-relativistic limit and the first relativistic correction for 
the nnvv process by expanding the propagator in terms of p/m 



± = _L( 1 + £f + o(4)). 

p.q moo \ mu x -m z/ J 
Application to the vector current amplitude yields 



M V = M V,NR 
V v 



AMj + 0(p 3 /m 3 ^ 



(22) 



(23) 



where the non-relativistic amplitudes M^ ,NR vanish and the leading corrections are given 

by 



AM 



GfCy 
2^2uj 



2 m 2 



fx (pi • q) - pi (pi • q) + {i <-> 2}^J T 



(24) 



A< = (25) 

with T nn the non-relativistic reduction of the 1 = 1 part of the T- matrix in Eq. (jll|) . The van- 
ishing of the non-relativistic vector amplitude generalizes the result of Friman and Maxwell 
[l], where this cancellation was observed for Landau-type interaction and OPE, to the com- 
plete nn T-matrix. This result is in fact, analogous to the absence of electric-dipole radiation 
in photon bremsstrahlung processes when the center-of-mass coincides with the center-of- 
charge of the radiating system, e.g. in pp bremsstrahlung. 
For the axial-current amplitude one obtains 

M, A = M*> NR + AM* + 0(p 2 /m 2 ), (26) 

where the non-relativistic amplitudes are given by 

M A ' NR = ^A[T nn ,S]; Mq' nr = 0, (27) 
2v 2u 

and the leading relativistic corrections are 

AM A = fr»* x (a ■ q) - BiT nn (p 1 ■ q) + {1 ±+ 2}) (28) 

AM ° A = ^Mt { Tnn(&i ' ^ " ■ Pi)Tnn + {l ^ 2} ) (29) 

with S = 8\ + S2 the total spin of the nn system. Eq. ()27|) has also been derived by Hanhart 
et al. Q|, and Timmermans et al. Q]. One sees from Eq. (|27|) that in the non-relativistic 
limit there is no contribution from the central interaction Tc, but the axial- vector current 
amplitude receives contributions from all other terms. 

The ppvv process can be treated analogously to nnvv process. The only differences 
are the coupling strength to the neutral weak current and the Coulomb corrections in the 
coefficients of the T-matrix. 

D. The npvV process 

In the n + p — > n+p + u + V process the momenta will be denoted by n and n' (p and p'), 
for the neutron (proton) in the initial and final state, respectively. In the ultra-soft region 
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(CO 



jp « 1) the vector current amplitude follows from Eq. (0) 



Ml 



G F 
2V2 



n.q n.q J v p.q p'.q 



J2 Fi np) \u(p')tt a u(p)u(n)n a u(n) + 



a=l 



'-yy <-> n'} 



and the axial- vector current amplitude from Eq. (jl(Jj) 



(30) 



Mi 



M A,atr + M 
v 1 v 



[A,dir 
v 

2mG F g A 



A,exch 
5 



2^/2 



La=l 



A+(n) 
2n.q 



-7,75^^fi Q )«(n) - {(n,nO - (p,j/)} 



+( 



u{n')Q a u{p)u{p')(Q a — 7^75 



2n.q 



+7^75 , P tt a )u(n) - {(n,n') (p,//)} 



(31) 



A+(p') 
2p'.g 

The exchange terms of axial- vector current matrix element are included explicitly in Eq. 
(I31|) . The direct part is analogous to the expression in Eq. (|21|) for the nnvv process. The 
only difference is the appearance of a minus sign in the term, where the neutron en proton 
momenta are interchanged. This is a consequence of the sign difference of the axial-vector 
coupling constants for neutrons and protons. 

The different structure of the exchange part (as compared to the nnvv process) comes 
from the sign difference between and c^. 

The expressions ()30|) and (j31|) simplify considerably, if one takes the non-relativistic limit. 
Using Eq. (J22j) one obtains for the vector current amplitude from Eq. (|3U|) 



M 



v 



a/2 2uj m ' uj 



and for the axial-vector amplitude from Eq. (|3*T|) 

Gf9a 



2V2u K 



(32) 



(33) 



where k = n' — n = p — p', D = {B\ — ^2), the spin exchange operator is P a = (1 + d[ ■ 0^2) /2, 
{..,..} denotes the anticommutator, T nP)dir and T np ' exch are given by the non-relativistic 
reduction of the direct and exchange parts of the np T-matrix. 




100 200 300 



p (MeV c" 1 ) 

FIG. 2: Cross section da/du; for n + n^n + n + u+ V as a function of neutron momentum in 
the cm. system, for u = 1 MeV, and summed over neutrino flavors. Shown are the result for the 
OPE (long-dashed) and the full T-matrix (full curve); in addition the separate contributions of the 
T-matrix Tti +Tt2 (short-dashed curve), Tso (dotted curve) and Tq (dashed double-dotted curve) 
are shown. 

Note that in order k/m « p/m in the np case there is a non- vanishing contribution for the 
vector current amplitude. This is analogous to the case of photon bremsstrahlung in NN 
scattering, where electric-dipole radiation is dominant for the np case. The commutator in 
Eq. receives contributions from tensor T T1 and T T2 , quadratic spin-orbit Tq and spin- 
orbit Tso components of the direct part of the np T-matrix. The anticommutator receives 
in addition to Tti, Tt2,Tq and Tso contributions from the exchange part of the np T-matrix 
also a central Tq contribution. 
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FIG. 3: Cross section da/dco for n + n^n + n + v + v as a function of neutron momentum 
in the cm. system, for u = 1 MeV, and summed over neutrino flavors. Shown are the result 
for the OPE(long-dashed), OPtRE (short-dashed curve), OPRSE(dotted curve), OPE without the 
"exchange" contribution (dashed double-dotted curve) and the full T-matrix (full curve). 

E. Comparison with one boson exchange(OBE) 

In this section we will calculate the neutrino emission cross section in free space. The 
expression for the cross section in the cm. system is 

,3 



da N f m i Jp 2 - luE + lu 2 /4 r „ , 

which is also given in Timmermans et al. |3j with the number of neutrino flavors Nf = 3. 
Neutrino pair bremsstrahlung has been calculated mostly, in Born approximation, with a 
two-nucleon NN interaction consisting of a long range one pion exchange (OPE) and a phe- 
nomenological Landau interaction as in Friman and Maxwell [l|. However, the use of lowest 
order OPE represents a severe approximation. First it is known that there is a substantial 
cancellation between the tensor contributions from rho and pion exchange. Secondly it is 
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FIG. 4: R = kiVi *- CTfll with 

a nr. the non-relativistc cross section and or the cross section, in which 
also the first order relativistic corrections are included. Taking the OPE as the NN interaction the 
nn and np bremsstrahlung ratios R are given by the dotted curve and the dashed curve, respectively. 
Also the nn bremsstrahlung ratio R is shown using the T-matrix for the NN interaction. 



questionable whether other (momentum dependent) interactions like the spin-orbit interac- 
tion Tso may be ignored. Hanhart et al. 2| found that the use of the full T-matrix leads to 
a reduction by a factor 4 compared to OPE for nn around saturation density. Our results 
for nn and np bremsstrahlung is a generalization of Friman and Maxwell's results: The 
amplitude is computed in terms of the (model independent) on-shell T matrix in stead of 
the Landau plus one-pion exchange interaction in the non- relativistic limit. The nn phase 
shifts are, for simplicity, assumed to be equal to the pp phase shifts, which are taken from 

In Fig. Elthe contribution of the various terms of the T-matrix T?\ + Tt2,T$o and Tq to 
the cross section in non-relativistic limit are shown separately for nn bremsstrahlung in free 
space. The contribution of the quadratic spin-orbit (the Tq term) and the tensor (the T?i 
and Tt2 terms) forces to the cross section cancel at low momenta. The tensor forces (the T?\ 
and Tt2 terms) dominate over the spin-orbit (the Tso term) and quadratic spin-orbit (the 
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Tq term) for momenta between 200 MeV/c and 300 Mev/c. From Fig. |2]one may conclude 
that at larger neutron momentum in the cm. system the spin-orbit force (the Tso term) 
becomes also important. 

Several results for the one boson exchange (OBE) contributions like one pion without "ex- 
change" contribution, pion, pion+tensor part of rho (OPtRE), pion+rho+sigma (OPRSE), 
exchange are shown in Fig. El as a comparison 29] . In the OBE potential contributions 
considered in this section the meson-nucleon form factors are not included. They can be 
neglected because of the relatively small momentum transfer, k < 2 p, involved. The OPE 
result overpredicts the full T-matrix result. At a neutron momentum of p ~ 300MeV/c 
in the cm. system the use of the the full T-matrix leads to a reduction of a factor of 4- 
5. Including the tensor part of the one rho exchange (ORE) to the OPE result is a much 
better idea. The cancellation of the tensor from OPE at short distance by the tensor from 
ORE, which has an opposite sign, leads to a result much closer to that obtained with the 
full T-matrix. The result for OPE without the "exchange" contribution, which is used in 
most "standard cooling scenarios" , is smaller than that for the full OPE, but has a different 
behavior than the result obtained with the full T-matrix. From a neutron momentum of 
250 MeV/c in cm. system the difference with the result of the full T-matrix increases. The 
contribution from one sigma exchange (OSE), which gives rise to a spin-orbit force, is also 
shown to give an estimate of the effect of the other mesons. The effect of the sigma is quite 
small. 

The calculations in Figs. El and El are done in the non-relativistic limit. Therefore it 
is important to check, whether the relativistic corrections are small. We can estimate the 
importance of the relativistic effects for OPE as well as for the the on-shell T-matrix taken 
as the iVW interaction. 

In Fig. 0]the relative relativistic correction R, in which the magnitude of the relativistic 
effects are compared to the non-relativistic cross section with OPE taken as iViV-interaction 
is shown. For the nnvv and the npvv processes the non-relativistic contribution comes from 
the axial- vector current. The relative relativistic correction R for the nnvv process remains 
below 15 percent and for the npvv process it remains even below 5 percent. Also in Fig. 0] 
R is shown for the nnvv process using the T matrix for the the NN interaction instead of 
OPE. The only contributions surviving the non-relativistic commutator in Eq. ()27j) come 
from the Ttx, Tp 2 , Tq and Tso parts of the T-matrix. The relative relativistic correction R 
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for T-matrix remains below 10 percent. Due to the chosen representation the spin-spin force 
is hidden in Tti,Tt2 and Tq. Some forces of the on-shell T-matrix have a non-relativistic 
character (scalar, spin-spin), while others don't have (tensor, spin-orbit, quadratic spin-orbit). 
In elastic scattering scalar and spin-spin forces dominate especially at low momenta. In 
the nnvv process these forces vanish in the non-relativistic limit of the bremsstrahlung 
amplitude, because they don't survive the commutator. They still have a non vanishing 
relativistic term in the bremsstrahlung amplitude, which explains the increasing importance 
of the relativistic corrections in the bremsstrahlung amplitude at very low momenta. 



III. NEUTRINO EMISSIVITY IN MEDIUM 



In this section we consider neutrino bremsstrahlung in a dense hadronic medium at finite 
temperature. In the simplest approach one can use the socalled convolution approximation 
(followed by Friman and Maxwell 1]) in which the free space bremsstrahlung process is 
folded with Fermi-Dirac single particle wave functions and the emission rate is obtained 
with the use of Fermi's golden rule. This approach is not applicable in more general cases, 
e.g. if one takes into account dressed propagators. To go beyond the convolution approach 
the more general framework of quantum transport theory Q, Q] is needed. The latter 
formalism and the application of the finite temperature Green functions is summarized in 
the appendix. 



A. The emissivity in quantum transport 

To compute the emissivity it is convenient to start from the Boltzmann equation (BE) 
for neutrinos (and anti-neutrinos), which schematically takes the form (see appendix A) 

d t + d q uj(q)d x 

where f v (q, x) is the single-time distribution function (Wigner function) of the neutrino with 
q the momentum and x space-time coordinate. The r.h.s. of Eq. (|35|) corresponds to the 
gain and loss collision integral (Appendix A and B). A similar equation holds for the anti- 
neutrinos. For a homogeneous system in Wigner representation the distribution functions 
become space independent. Furthermore the time dependence of the collision integrals can 



f u (q, x) = I v + (q, x) - J+ (q, x), (35) 
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be neglected. Therefore we drop the x argument at the r.h.s. of Eq. (|35p. The use of 
the BE provides a general formalism for neutrino and anti-neutrino emission, absorption 
and scattering. The collision integrals I h and I + ~ are directly related to the neutrino 
selfenergies $ h and $ + ~ (Eq. (|B2p ) which in turn are expressed in terms of the hadronic 
polarization S~J~ ,+ ~(q) and the leptonic couplings and propagators (Eq. (|B1|0 . The former 
are closely related to retarded polarization or the current-current correlation functions 

S?{q) = SjT(-ff) = 2ig B (u;)^mU^(q) = Am J d 4 £exp (zg£)<^(0«£)> (36) 

with the retarded polarization function U R (q). The general polarization receives contribu- 
tions from vector, axial-vector and interference terms 

n^( ff ) = c 2 v ulM + 4n^( g ) + c^n^fe), (37) 

where H^fa), H^ v (q) and UY^(q) are the vector, the axial-vector and the mixed part. In 
general one has four independent components Il^ (g),n^(g),n A (g) and U VA (q) |17j . The 
lepton couplings and propagators give the leptonic tensor = 8(q± q 2 + q1q 2 — { < lx- ( h)9 tl ' ,v — 
ie a ^q lM2>a ). 

In the present case of emission we take the neutrinos to be free. The emissivity (the power 
of the energy radiated per volume unit) is obtained by multiplying the energy with the l.h.s. 
of Boltzmann Equation (BE) (see Appendix) for neutrinos and anti-neutrinos, respectively, 
summing the neutrino and anti-neutrino expression, and integrating over a phase space 
element: 

d f d 3 q 
euP = JtJ (2<j 

From Eq. (|33|) follows 

d 3 q 

where I^ + ' cm (q) and l£~' em (q) are the terms of the collision integrals, which correspond to 
neutrino emission process. 

To obtain the emissivity the leptonic tensor has to be contracted with the structure 
function 

e„p = -2X; 

/ 

(27r) 4 5 3 (<fi + q 2 - q)8{uj{qi) + u{q 2 ) - u) [uj{q\) + uj(q 2 )) 



j^[fMt) + fMt)Mq)- (38) 



f d?q 2 


f d 3 q! 


f d 4 q 


1 (2n) 3 2io(q 2 ) J 


(27r) 3 2cj(gi) J 


(27T) 4 



g B (u)A^( qi ,q 2 )ImU R (q). (40) 
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The number of neutrino flavors is included by the summation over f. For neutrino pair 
bremsstrahlung it is more convenient to use in the leptonic tensor q = q± + q 2 instead of q± 
and qi. Using Lorentz covariance, we can write 
d 3 qi d 3 q 2 



J CJi LUo 



= 5 (?V - <zV") I S\q -q 1 -q 2 )= 1 -^ (gV - <?<T) • (41) 

This simplifies the expression of the emissivity 

^^E/^-WM (42) 
with W(q) = -2 g B (uo) L^(q) Smll* (g). 



B. Hadronic polarization 

Which type of correlation diagrams are dominant in the neutrino-hadron interaction 
processes depends strongly on the kinematics. In particular in the space-like region (|g| > 
u) (scattering) the one-loop QPA diagram and its random phase approximation(RPA)- 
type iteration dominate; in contrast in the time-like regime (u > \q\) the QPA process 
is kinematically forbidden and two-body (and many-body) collisions are required as was 
already clear from the discussion of the free space case. 

For practical calculations of the polarization one has to make a choice between the use of 
dressed Green functions and the use of quasi-particle Green functions. On the one hand the 
use of QPA in the soft limit of uj —>■ 0, leads to the property that Smll^ behaves as 1/u 2 in 
all orders, i.e. an infra-red divergence (this behavior is correct only for the free case, where 
the external legs are on-shell). Hence one expects that in the soft limit non-perturbative 
effects play a role (see the LPM effect, below). On the other hand as pointed out in ref. 
in using dressed propagators special care has to be taken to avoid double counting, i.e. one 
has to restrict oneself to socalled proper "skeleton diagrams" . (An example is the two-loop 
self-energy insertion in diagram which is already effectively included in the one-loop 
diagram with full Green functions.) Another problem connected with the use of dressed 
Green functions and bare vertices is the conservation of the vector current. 

In general an expansion in terms of QPA diagrams is simpler (than in terms of full 
Green functions) since there are no such spurious diagrams, and also current conservation 
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+ 




FIG. 5: One loop contribution to S+, 

is satisfied at each loop level. Below we will show that in the special case of a imaginary 
part of the self-energy (width) there exists a 1-1 correspondence between the QPA and 
dressed Green functions diagram expansion, i.e. the proper full diagrams can be expressed 
as multiplicative correction factor u 2 /(u 2 + T 2 ) to the QPA result. This result allows us to 
use the QPA and include the finite width at the end. 

At low temperatures the leading diagrams are those which contain a minimum number of 
off-diagonal G + ~ and G h . In the closed diagrams the +- and -+ lines are cut. In the QPA 
limit this gives back the original Feynman graphs. The + part is the Feynman amplitude 
and - part belongs to the conjugated Feynman amplitude. 



1. One-loop in QPA 

For completeness we give the one loop polarization function in the QPA limit 

iS?{q) = / ^^jAT»G + (p)T u G^(p^m q + p' -p), (43) 



where T M = fjk-y^cv - c A js) 



In the non-relativistic QPA limit Eq.(j43J) can be factorized in terms of a hadronic loop 
and couplings X^ v 

d 3 p d 3 p' 
(2tt) 3 (2tt) 3 ' 

(2tt) V(g +p'- p)X fJlU = 2g B (u)X f ,J (q) (44) 

with 

^ 

After integration I (q) = f^£(<?) witn |l3 

C(q) = ln((l + exp{-/3(e_(g) - ^)})/ ln(l + exp{-/3(e + (g) - /i)})), (46) 
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FIG. 6: The 3 different types of "closed diagrams" at the two loop level. These diagrams can 
be considered as a) (lowest order) propagator, b) vertex and c)interaction renormalization of the 
QPA. 

where e±(q) = (lu 2 + e|)/4e 9 -± lu/2 with e 9 = <f /(2m*). 

One sees that in the one-loop approximation in the QPA only the space-like contribution 
(uj > \q\) does not vanish. 



2. Two-loops in QPA 

In Fig. injtlie 3 different types of "closed diagrams" at the two loop level are shown. These 
diagrams can be considered as (lowest order) propagator, vertex and interaction renormal- 
ization of the one loop in QPA, respectively. We begin considering the simple case of nn 
neutrino pair bremsstrahlung with the on-shell T- matrix in Eq. ((UJ. Diagram contains 
terms with a causal propagator G ++ and an acausal propagator G with the same argu- 
ments,which can be pi — q or p\ + q, whereas diagram E)d contains terms with a G ++ and 
G with different arguments (opposite signs for q) and one obtains for diagram EK + Eb 

5 5 - - 2 d A Pi d% i d^k 



i(s^\ q ) + = EE/[I1 

1 n i « „■ 1 



a=l 0=1 ' 



lrpl* 



\ (2tt) 4 (2tt) 4 J (2tt) 4 a 13 
(Tt[Q-^G + ( P2 )Q^G^(p' 2 )] 

Tr[A-^ 1 G + (p l )A+X 1 Gt(p l )] + {1 - 2}] 
(2vr) 8 5 4 (A; + p' 2 - p 2 )5\q + p[ - k - p x ) (47) 

with = ^ + G ++ fe-g)r, + r,G ++ (^ + g)fi+ + ,A~ = aud V, = ^(c v - 

C475). The definition of Q ++ is given in Eq. ()12j) and Q follows from the relation Q = 
7 f2 ++ 7 . Thus in the non-relativistic limit diagrams EK andEb have opposite signs ±l/u>. 
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One obtains for diagram |Ht 

5 5 - 2 d 4 Pi d A Pi i d 4 k 



^ + ' (c) = /ee[ii 



Irpl* 



rTzT, 



4 ± a ± (3 

Tr[A-- 2 G^(p 2 )n+ + Gr(P2) 
TV[0--^G-+( Pl )A+X 1 G+-(p' 1 )] + {1 - 2}; 
(2vr) 8 5 4 (A; + p' 2 - p 2 )<5 4 (g +p' 1 -Ar- Pl ). (48) 

Note that only the -+ and +- lines are cut in the diagrams of Fig. |3 since cutting the 
T-matrix would lead to double counting. 

The above expressions become simpler, if the QPA Green functions are used (see Eqs. 

flsani ). 

(2tt) 8 5 4 (A; + p 2 - p 2 )5 4 (g + ^ - k - p 1 )X fiv , (49) 

where X contains all operators and f(Ei) = (exp ((3(Ei — //)) + l) with the relativistic 
energy and /i the relativistic chemical potential. In particular for diagrams IBH and |Hb we 
obtain 

x<£ + x® = EE^^Trr^^n-^^)] 

a=l (3=1 

Tr[A+; i A+(p*)A^- 1 A + (p'*)] + {1 <-> 2}), (50) 

and for diagram Efc 

= EE^T;(Tr[A5' a A + ( P ;)(]-V(p';)] 

a=l /3=1 



Tr[fi+ + A + (pDA,;,i A+ (^)] + i 1 ^ 2})- ( 51 ) 

One verifies that the sum of all two-loop diagrams conserves the vector current,i.e. 
q^S~„(q) = 0. First diagram IHfc is current conserving on its own. That the sum of dia- 
grams ink and IHb is current conserving can easily be deduced from Eq. (J2U)| by noting that 
Q^iPn/ (Pi ' q) ~PiJ (Pi ' <?)) = 0- 111 the following the hadronic part of interaction matrix X^ v 
is evaluated in the non-relativistic limit for the cases A, V, V A separately. 
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a. The non-relativistic limit Although in principle X can be evaluated relativistically, 
we will use the simpler non-relativistic formalism. First we consider the vector current XY V . 
Expanding the Green functions G ++ and G (see appendix IA"|) in powers of (p- q)/{m*uj) 
leads to 

X ( £' V + X$' V + = C ^V,V U \T™\ 2 + 0(\p\ 3 /m*% (52) 



where 



and 



V, = -A- ( - + + V'i»{l + ^) + {1 - 2}) (53) 



\T nn \ 2 = 4(\Tc\ 2 + \T Q \ 2 + \T Tl \ 2 + \T T2 \ 2 + 2\T so \ 2 ). (54) 

We see that in leading order in the non-relativistic limit, with G ++ (p±q)r^ — > ±p^/(m*uj), 
the vector contributions cancel due to p' x +P2^ —p'ia ~V2 i i ~ 0; while the separate diagrams 
do not vanish. 

For the axial-vector current to obtain the non-relativistic limit we expand the G and 
G ++ functions in powers of (p-q)/ (m*u), and replace the coupling —> a-p/m* S^o + a^^i. 
For diagrams + Eb the hadronic part of the interaction matrix is 

X^ A + x!S' A ^0(\p\ 2 /m* 2 ), (55) 

x^ a +x^ a = ^j:m 2 {^[{ti*tm><i) j \ 

+ Om 2 /m* 2 ))+{l^2}) (56) 



(Pi -P'l) ■ (1 , n ,n2 / 2 



and 



rrru) 



+0(|p] 2 /™* 2 ) + {l «2}) (57) 



with 



^ = (Ji, ^2, *4, ^) = (ft, n, k, k', ft). (58) 

Hence in leading order in the non-relativistic limit there is a contribution to the axial- vector 
current (the contribution from the central interactions to the axial-vector current vanishes, 
since they commute with the weak spin operator). For diagram U one has 

X^ A ^0(\p\ 2 /m* 2 ), (59) 
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x ^,a = C -^fJ2J2l6(^Tr + T u %^ 

*U v =2 u=2 

ft x i u )% x i v y(i + ^+^+^+^)-g + 0(bf/m ^ (60) 



^ j;=2tt=2 



ft x Z^fe ■ ((pi+£ +^ +^ 2 ) x 4)) + 0(\p\ 2 /m* 2 ) (61) 

with z = j = 1, 2, 3. The definitions for T are given in Eqs. (|18|). 

As for the mixed VA contribution in the non-relativistic limit the traces vanish and hence 

x (a),VA + X (6),VA + x (c),VA ^ (p 2 /m* 2 ). (62 ) 

Therefore we obtain the (well known) result |l( that in leading order in the non-relativistic 
limit there is only a non- vanishing contribution from the axial vector current. In free space 
Fig. |3 shows that at momenta relevant at nuclear matter densities the p/m corrections to 
neutrino pair emission are only of the order of 10 percent. Therefore one may conclude 
that the leading order non-relativistic result with only the axial-vector current constitutes a 
good approximation. Contracting the polarization function S~J~(q) with the leptonic tensor 
L^(q) = if (gV - gV A ) yields 



W(q) = i^[L^{q){S-^\q) + S^ b \q) + S^\q))] 
_ 2ng 2 A G 2 F f 2 r d 3 Pl d 3 p\ 

\M\ 2 5 A (q + p' 1 +p' 2 -p 1 -p 2 ), (63) 



where 



\M\ 2 = 32 ]T |T;| 2 ((2c 2 - \q\ 2 )\l v \ 2 - (q- l v , 

v=2 

2/, .2 



+!6 E E (W + TvX 1 ) ( (ft x 4) ■ qf + ft x f w 

v=2 u=2 \ 



)V-ir) • (64) 



The LPM-effect 



The free space AiV neutrino-pair bremsstrahlung process exhibits an infrared 1/u di- 
vergence (see section HTC} . The QPA result in Eq. (jB3J) also shows an infrared divergence, 
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reminiscent of the free space bremsstrahlung. It is well known that the singularity in the elec- 
tromagnetic bremsstrahlung process is suppressed in a medium,the Landau-Pomeranchuk- 
Migdal(LPM)-effect, whenever the mean free path of the emitting particle becomes compa- 
rable to the photon formation length, u — v-q. The former is characterized by the imaginary 
part of the self energy, r, of the emitting particle, while the photon formation length can 
be approximated in the non-relativistic limit by the formation energy, to. Therefore the 
LPM-effect is expected to become effective whenever u w V. 

The LPM effect has been discussed recently in various contexts. For instance for photon 
emission in a quark gluon plasma by Aurenche et al. ^| and Cleymans et al. in terms 
of thermal field theory. Analogously one expects similar effects in the electroweak case (as 
noted by Raffelt (^(J for neutrino pair and axion production in supernovas and neutron 
stars). 

Here we estimate the LPM effect on the response function S and the emissivity as a 
function of the temperature and density by using the dressed propagators for G~ + ,G + ~ and 
G in Eqs. (jCllllC13j) . Note that in a fully dressed Green functions formalism diagram a) 
of Fig. |B1 is not proper skeleton diagram and its contribution is already included in the fully 
dressed one-loop diagram. In this case the appropriate irreducible diagrams to be consid- 
ered are given by the dressed one-loop, the corresponding two-loop vertex correction (these 
together conserve already the vector current) and the two loop interaction normalization. 
We evaluate these in the limit of Y — 23mS < 9fte£. Following we note that in the limit 
T = constant and q — > it is possible to write the fully dressed diagrams in terms of the 
lowest non- vanishing order in the QPA in the low temperature limit. A remark has to be 
made about the one loop result 

We note that it is possible to relate the one loop result to the lowest non- vanishing order in 
the QPA, if the quasi particle width T in the numerator is represented by the one loop QPA 
self energy 



3m £ Q 



(66) 
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In the low temperature limit one can make the approximations f(E p > + to) ~ and f(E p — 
which leads to 




We note that with the present (which differs from the result given in p|, namely 
= T^fef^ ) the CVC relation holds, because the vector current is conserved in the 



QPA limit. 

If only the dressed off-shell propagators, G ++ and G are kept while G + ~ is replaced 
by Gq~ we find a different result for the damping, C = jt^tji- From this we conclude that 
the dressing of all G"s should be considered on equal footing. In some previous works (e.g. 
Raffelt and Seckel 20]) the quasi-particle width has been included directly (in a rather 
adhoc fashion) in the cross section by replacing, by a modified one, ^-^cPr 2 ' wriere a is 
taken to be unity. In this case V is purely a parameter with no microscopic origin; in reality 
T depends on momentum, density and temperature. 



D. Modification of the T-matrix in the medium 

Above we have considered the T-matrix in free space. In the past the possible medium 
modification of the T-matrix has been addressed only in a very few papers. The Rostock 
group has studied the effect of the medium on neutrino emissivities |22| in the frame work of 
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a thermal dynamic T-matrix. It was found that at T = 4 MeV the ratio R of emissivities for 
the in-medium T-matrix to the free T-matrix result is about 0.8 for nuclear saturation density 
for the modified URCA process, and a striking m 0.05 for the neutral current bremsstrahlung 
process. The latter effect was ascribed to the Pauli blocking of the low momentum states. 
The results were obtained using a separable approximation to the potential neglecting 3 P 2 — 
3 F 2 tensor coupling. Here we estimate the medium effect by using a G-matrix at zero 
temperature to account for Pauli blocking, which includes the full tensor force. The Bethe- 
Goldstone equation for the G-matrix is 

G<jt,0> = V(f,f) + £ / gfrVV (70) 



A,/ 



where X,i are the helicities and isospin of the intermediate state, respectively. The single 
nucleon energy above and below the Fermi momentum pp are E(p) and e(p), respectively. 
Here the G-matrix of Banerjee and Tjon [3] in the lowest order Brueckner theory (LOBT) 
is used; the single nucleon energies are given by 

E{p) = f- ; e(p) = A + ,r 



2m 2m* 

The gap A and the effective mass m* are determined in the LOBT in a self consistent way. 
As an interaction in Eq.( [7UJ) the Bonn-C potential is used and for Qp au ii m Eq.( 170)1 an 
angle averaged Pauli operator is used to construct the G-matrix. 



IV. RESULTS AND DISCUSSION 



We will compare the emissivity of the nn neutrino pair bremsstrahlung for the different 
NN interactions at densities n = l/2riQ,no and 2no at T = lO 9 ^. We will derive the 
expression of the emissivity in the non-relativistic QPA limit starting from Eq. (|42|). From 
Eqs. (I63j) and (|64j) the function W(q) is obtained. In here the momentum q is neglected 
in the momentum conserving delta function, because it is much smaller than the neutron 
momenta. Next we separate the angular and energy parts of the nucleon phase space by 
performing the angular integrals with the momenta of the degenerate neutrons approximated 
by the neutron Fermi momenta. Finally we use the independence of the matrix elements of 
q and introduce the dimensionless parameter y = uo/T to simplify the expression and one 
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neutron matter 


symmetric matter 


density 


1/2 n 


riQ 


2 no 


1/2 n 




2 no 


m* jm 


0.77 


0.64 


0.49 


0.66 


0.58 


0.49 


OPE 


7.3 


4.8 


2.1 


2.7 


2.4 


1.6 


OPRE 


3.9 


2.3 


1.0 


1.7 


1.3 


0.8 


OPROE 


3.2 


2.0 


0.9 


1.4 


1.0 


0.7 


OBE 


3.7 


2.5 


1.3 


1.5 


1.2 


0.9 


OPE+TPE 


10.2 


6.9 


3.7 


3.5 


3.2 


2.3 


OPRE+TPRE 


4.4 


2.8 


1.4 


1.8 


1.4 


1.0 


OPROE+TPROE 


1.2 


1.2 


2.7 


0.4 


0.3 


0.4 


OBE+TBE 


1.4 


0.6 


0.3 


0.8 


0.5 


0.2 


T-matrix 


2.2 


1.1 




1.1 


0.7 


0.4 


.R-matrix 


2.5 


1.4 




1.2 


0.8 


0.5 


G-matrix 


2.7 


1.6 


0.6 


1.2 


0.9 


0.5 



TABLE I: emissivity in 10 erg cm s 



obtains 



where 



4G 2 F g 2 A m**p Fnrr8 



15(2tt) 



T 8 J dyd cos(6 , i 2 )<icos( 



H{s,t) 



V2 + 2cos0i2 



J(v) 



(47rV + n 7 ) 
6(1 + exp(y)) 

and the hadronic part of the interaction matrix 

H( s ,t) = ( 8 j:\r v ^ + 2Y: £ (r^ + xrrj 

u=l «=2 v^u;u=2 



(71) 



(72) 



(73) 



is a function of the Mandelstam variables s, t and pi? n is the neutron Fermi momentum. The 
integration variable 12 is the angle between pi and p 2 and Q\y is the angle between pi and 
We note that in the limit that the G-matrix is replaced by the anti-symmetrized one-pion 
(or one-rho) exchange potential Eq. (|71|) reduces to the result of Friman and Maxwell [l|. 

The results for the emissivities are summarized in table|l]for neutron matter for 3 different 
densities, n = l/2n ,n and 2n at T = 10 9 K. It is seen that (similarly as in free space) 
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compared to the T-matrix result the anti-synimetrized OPE overestimates the emission rate 

n 

by roughly a factor 4; this is in agreement with the conclusion by Hanhart et al. [2|. If the 
exchange terms in OPE are (arbitrarily) omitted the result is close to that of the T-matrix. 
In the past in some cases phenomelogical correction factors are also introduced to simulate 
initial and final state interactions as a correction to OPE , which tend to reduce the 

OPE result. Contrary to naive expectations based on upon the Pauli blocking mechanism 
we find a slight increase in the rate if the T-matrix is replaced by the in-medium G-matrix 
as calculated by the Bethe-Goldstone equation described in the previous section. 

To obtain more insight in the medium effects we also listed in the table the separate results 
for one pion, pion+ rho, pion+rho +omega (OPROE), and full OBE (which includes also 
sigma exchange), and also for the corresponding OBE plus iterated OBE (i.e OBE+ TBE). 
It is seen that one rho-exchange gives a substantial cancellation of the OPE (also observed by 
Friman and Maxwell 1]); on the other hand the iterated OPE (referred to as TPE) leads to a 
stronger tensor force, and hence a larger rate. It is also seen that the contributions of omega 
exchange and sigma exchange (which contributes mainly to the spin-orbit interactions) are 
non negligible in particular in the TBE process, and as a consequence in the combined 
OBE and TBE contribution becomes even smaller than the full T-matrix result. These 
momentum dependent interactions do not appear in the conventional Landau fermi liquid 
interaction, but do seem to play role in the weak bremsstrahlung. We attribute the finding 
that the G-matrix gives a slightly larger contribution than the free T-matrix mainly to a 
Pauli blocking of the TBE contributions, and hence a smaller destructive interference. We 
note that the present result deviates from the one in Q where at „ = n in neutron matter 
the ratio of the rates computed with G-matrix and T-matrix was found to be 0.05; a possible 
explanation could be the neglect of the 3 P2 — 3 i<2 tensor coupling in that work. 

As to the density dependence the decrease of the rate with increasing density is mainly 
caused by the variation of m* in Eq. (fTT|) and to a lesser extent by the different ranges of 
the various meson exchanges. For completeness we also show the corresponding results for 
symmetric nuclear matter in table H] The weaker density dependence in this case can be 
attributed to less variation of m*. 

Finally we turn to the LPM effect. Clearly its possible relevance in the present case 
depends on the magnitude of the width T, which is a function of u, temperature T and 
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y 



FIG. 7: The functions f(y,T) and Ilpm(u,T) are shown at temperature T = 2 MeV (dotted 
curve), T = 5 Mev (dashed curve) and T = 10 MeV (dashed-dotted curve). The QPA result is 
given by the solid line. 
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density. Here we use the parameterization 



r(^,T) = a(^ + T 2 ) (74) 
to be able to estimate the importance of the LPM-effect. For u < 60MeV and T < 20MeV, 



this roughly coincide within a factor 2 with Aim et al. 



Including the LPM effect in the 



emissivity the function I(y) in Eq. (ffT| has to be replaced by 

I L PM(y,T) = y 2 f(y,T)I(y) (75) 

with f(y,T) = l/(y 2 + r (^ r ) ) ; which can be derived from Eqs. (jHTjl . (}68|) and (}69j) . The 
function f(y,T) describes very roughly the behavior of Smll^. To give an indication of 
the importance of the LPM effect and to demonstrate the influence of the weighting factor 
in the emissivity we show in Fig. [7|how the functions f(y,T) and lLPM(y,T) in Eq. (|75Jl 
are modified for various values of the temperature. The value of the parameter a depends 
weakly on the density and is approximately 0.2MeV 'c~ x . One sees that the function f(y, T) 
has a singularity at y — in the QPA. The LPM effect suppresses this infrared divergence. 
The function IiPM{y, T) is less sensitive to the LPM effect compared to the function f(y, T), 
because the weighting factor in the emissivity strongly suppresses the y = contribution. 
Therefore the LPM effect in the emissivity is negligible for T < 5 MeV. Comparing the ratio 
of the emissivity with and without LPM effect Rlpm = t/tipm at T = 5 MeV,T = 10 MeV 
and T = 20 MeV gives 0.89,0.68 and 0.35, respectively. The influence of the LPM-effect in- 
creases with temperature and becomes appreciable above T = 5 MeV. Therefore in practice 
in calculating the emissivity the LPM effect does not play an important role for small T, 
say T < 5 MeV. 

Finally we note that an additional medium effect, not considered here, is the possible 
medium effect of the axial vector coupling g^, which has been considered in j^j], where it 
is was found that space like axial coupling is quenched by about 20 percent. However the 
timelike axial coupling is not necessarily equal, since Lorentz invariance is broken. Experi- 
ments with first-forbidden /3 decay of light nuclei give indications for an enhancement of the 
time-like axial charge of about 25 percent in the medium |27j. This is in agreement with 
meson exchange calculations in the soft pion approximations 
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V. SUMMARY AND CONCLUSION 



In this paper we studied the neutrino emissivity for the neutral current NN 
bremsstrahlung process, relevant for neutron star cooling. In particular we considered some 
effects that are not included in the standard cooling scenario of [1], which is based upon a 
non-relativistic quasi-particle approximation and the use of the one-pion exchange potential. 
The effects considered, namely the description of the NN interaction, the LPM effect and 
relativistic effects, influence the neutrino emission of the neutral current bremsstrahlung 
process. Therefore these effects are expected to affect also other neutrino emission processes 
in a similar way. 

First we studied how the description of the NN interaction influences the NN 
bremsstrahlung process. In the low density limit using the fact that to is small the Low 
theorem \^ can be applied, which allows us to use the on-shell T-matrix, specified by em- 
pirical phase shifts, and to compare it with OPE. At typical neutron momenta in neutron 
stars, approximately 300 MeV/c, the resulting free space cross section is roughly a factor 4-5 
reduced compared to the application of OPE. Although adding ORE to OPE is an improve- 
ment, the result still differs a factor 2-3 with that obtained using the T-matrix. We also 
analyzed which Fermi components of the T-matrix dominate the rate, namely the tensor 
and the spin-orbit type terms. 

To evaluate neutrino-pair bremsstrahlung in a finite medium at finite temperature we have 
used a closed diagram technique up to two loops. It is found that at n ~ n the neutrino 
emissivity, applying the on-shell T-matrix to describe the NN interaction, is roughly a factor 
4 smaller than those based upon OPE. This is in qualitative agreement with the conclusion 
of Hanhart [2] . Including medium effects from Pauli-blocking by replacing the T-matrix by 
a in-medium G- matrix, we find a small increase of the emissivity of 20-30 percent. 

Secondly in order to investigate the many-body correlations we can go beyond QPA by 
considering dressed propagators with a temperature dependent imaginary part T. Of course 
gauge invariance of the vector current is conserved in our approach. In particular we find 
that in the medium the damping of the infrared divergence, the LPM effect, has a negligible 
effect for low temperatures (T < 5 MeV) ; this is due to both the small single-particle width 
(r « T 2 ) and a weighting factor depending on uj in the phase space integral. 

Finally we estimated relativistic (recoil) effects to be rather small, of the order of 10 
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percent, at nuclear saturation densities. 

In short the description of the NN interaction by the on-shell T-matrix OPE has the 
largest impact on the neutrino emission of the bremsstrahlung process; roughly a reduction 
factor of 4. Other effects are relatively small; below 30 percent for T < 5MeV. 
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APPENDIX A: NEUTRINO TRANSPORT 



In the present paper we use the finite temperature real time Schwinger-Keldysh formalism 
to compute the collision integrals in the transport formalism. For the sake of completeness 
the main steps are summarized in this appendix; for more details we refer to ffj. In this 
formalism one must distinguish between vertices with indices (+) and (-). For given real 
interaction these are associated with the value —iV (time ordered part) and with adjoint 
vertex +iV (anti-time ordered part). The corresponding finite temperature Green functions 
(applied to neutrinos as well as the nucleons) can be expressed as a two times two matrix 
propagator: 

Sometimes it is more convenient to use the retarded and advanced functions: 

iG* 2 = e{t x -* 2 K{vK*i),fe)}), iGf 2 = e{h - t 2 K{v^i),fe)}>, (A2) 

The propagators satisfy the Dyson equation 

G(xi, x 2 ) = G (xi, x 2 ) + Goix!, x 3 )<t(x 3 , x 2 )G(x 2 , Xl ) (A3) 

where $ is the proper self-energy. 
Equivalently in integro-differential form 

@iQ lj2 = S 1;2 a z + a z J rf3$ 13 G 3i2 (A4) 
$ 2 G lt2 = 8 x , 2 a z + J d3G li3 ^ 2 a z , (A5) 

where g_ z is the Pauli spin matrix. The semi-classical neutrino transport equation are ob- 
tained by subtracting the Dyson Eqs. (|A4|) and (|A5J) for d\ and d 2 

iG{x u x 2 ) ft X2 -i $ Xl G{x u x 2 ) 
= G(xi, x 3 )®(x 3 , x 2 )a^ - OjMxu x 3 )G(x 3 , x 2 ), (A6) 

In particular the transport equation for the off-diagonal matrix Green function reads 

-Me<S> R (x 1 ,x 3 ),G + -'-+(x 3 ,x 2 )} - [toeG R (x u x 2 ),$ + ->- + (x 3 ,x 2 \ 
= \ {G + ~'- + (x h x 3 ), $ + -'- + (x 3 , x 2 )} + l - {$ + -'- + (x 1; x 2 ), G + -'- + (x 3 , x 2 )} , (A7) 
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As a result of the assumption of the existence of the Lehmann representation we have 
UeG R = %teG A = KeG and Ke$ K = Ke$ A = Ke$. The Wigner transforms of the off- 
diagonal Green functions correspond to Wigner densities in 4-coordinate and 4-momentum 
space. In the gradient expansion the Wigner transforms of convolution integrals can be 
expressed in terms of Poisson brackets {A, B} PB = dkAd x B — d x AdkB. This leads to the 
quasi-classical neutrino transport equation in which the neutrino self-energies enter in the 
loss and gain terms 

i {^e G~\p, x), G + -<- + (p, x)} pB + i {3?e G(p, x), (p, x)} pB 

= G + ->- + (p, x)$ + ->- + (p, x) + <5> + ~>~ + (p, x)G + ->- + (p, x), (A8) 

The first Poisson bracket at the l.h.s. side leads (Vlasov part) to the Boltzmann drift term, 
whereas the second one corresponds to off-mass shell effects. After separating the pole and 
non-pole terms: 

G+->- + (p,x) = G+-'- + {p,x) + Gt s >- + (p,x) 
the quasiparticle part of the transport equation is given by 

i{^eG-\p,x),Gt->- + {p,x)} pB = G + -<-+{p,x)<$> + -<- + {p,x) 

- $ + ->- + (p,x)G + -'- + (p,x), (A9) 

where 3ieG~ 1 (p, x) = ftp — $le& R (p,x). The l.h.s. corresponds to the drift term of the 
Boltzmann equation and the r.h.s. to the collision integrals. The remainder part of the 
transport equation 

i{^eG- 1 (p,x).Gtf f '- + (p,x)} pB + l {^eG(p,x),^- + (p,x)} pB = (A10) 
describes the off-shell effects, which we neglect. 

The on-mass-shell neutrino propagator is related to the single-time distribution functions 
(Wigner functions) of neutrinos and anti- neutrinos, f u (q) and fp(q), 

G + (q, x) = ^ [5 (g - oo(q)) f v (q, x) - 5 (q + u(q)) (1 - f p (-q, x)) } (All) 
uj{q) 

and for the Gq~ propagator f v (q) replaced by 1- f u {q) and 1 — fp(q) by 1 — fp(q). In this 
limit the Boltzmann equation for the neutrino distributions is obtained 

d t + d q u(q)d x ]f u (q,x) = J™^Tr[<S>- + (q,x)G+-(q,x) 

-^ + -(q,x)Go + (q,x)]=l; + (q,x)-ir(q,x), (A12) 
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where the r.h.s. corresponds to the gain and loss term, (the Boltzmann Eq. for anti-neutrino 
follows by integration over the negative go)- 

APPENDIX B: COLLISION INTEGRALS 

In the lowest (second) order in the weak interaction the neutrino transport self-energies 
are given by 

- (g,x) = J ^^(2vr)V( gi + q 2 - q)iV^ zG» + (q 2 , x)zT^S^ + - (q u a#31) 

where S~\' ~(q) 1S the baryon polarization tensor, and is the weak leptonic interaction 
vertex. 

The collision integrals in Eq. (jA12j) . which are expressed as a convolution of the lepton 
self-energies $ and the intermediate (anti-)neutrino propagator, consist of a sum of a loss 
and a gain term; e.g. the neutrino gain part 

I- + (g» = r^Tr[$-+(g,x)G + -(g,x)] (B2) 
JO Ztc 

contains a (space-like) scattering (proportional to f v (l — f u )) and a (time-like) pair emission 
term (oc (1 — /p)(l — f u )) The anti-neutrino one is obtained by replacing the positive energy 
range by the negative one. 



APPENDIX C: FINITE T HADRONIC GREEN FUNCTIONS 

Although in the neutrino sector the stationary condition $ h G+~ = $ + ~G~ + is not 
satisfied (see appendix), in the hadronic sector it is. Therefor the nucleons can be treated in 
the equilibrium Green function's formalism. The retarded self energy E R can be decomposed 
in Lorentz components, in nuclear matter only the scalar and vector components are non 
zero 

E*(p) = £f(p)+*#(p) 
with p = (p°,p). The retarded relativistic dressed baryon Green function [l^ is 
r R (n s = ^ + m-^(p) + Sg(p) 
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and the spectral function 

A{p) = -2%mG R {p). (C2) 
Using Eqs. ICJ1I and IC2I we can now give the following relations 

G- + ( P ) = if(p°)A(p), (C3) 

G + -(p) = -z(l-f(p°))A(p), (C4) 

G~(p) = (1 - f( P °))G R (p) + f( P °)G A (p), (C5) 

G ++ (p) = -(1 - - /(p°)G B (p) (C6) 

with /(po) = l/( ex P(/5(P° — /-0) + l)j = and the chemical potential /x = + 
3te££ w (p F ). We will now define the relativistic effective Dirac mass mp = m + 9fte£f (p), 

pO = p o - KeE^Go), p = p + 3?e£$(p) ,E P = ^{p) 2 + {m D ) 2 and T = 23m( - - 



^T,g(p) + We will consider two cases:i) the QPA Green functions: SmS(p) — ► 

and ii) the non-relativistic Green functions. 



1. Green functions in QPA 

In the QPA case, the imaginary part of self energy 9mS(p) vanishes. This gives the 
following definitions for the Green functions in Eqs. (|C3Jl . (|C4Jl . (|C5J) and (|C6jl 

Gt-(P) = -2m^M(l - f(p°))5(p° - E p ), (C7) 
G + (p) = 2z7T^|^/(p )5(p - E p ), (C8) 

G -(p) = (G ++ (p)r = ^^P, (C9) 

where we have the positive-energy operator A + (p) = and Zp X (p) = d(p 2 — m 2 D )/dp°. 
The causal propagators Gq ~ and Gq + are off-mass- shall. If p is on-mass-shell p 2 — mo 2 = 0, 
then G can be rewritten as 

Grip ± ff) = (G ++ (P ± (?))• = 2m * A _ + ® . (CIO) 

We point out that, when taking complex conjugates, it is understood that Dirac gamma 
matrices are not conjugated. The free case can easily be obtained from this. By replacing 
mo,p by m and p we obtain the free Green functions. 

34 



2. The non-relativistic Green functions 



In this part will be given the non-relativistic Green functions. Besides the non-relativistic 
limit, we will assume that the width of the quasi-particle state, is small, SmS fl (p) <C 
5te£ fi (p). We will now define the Green functions in the non-relativistic limit as 

G+ ~ M= w^wTwm (1 - m) (C12) 

and 

G -- M . (G++wr . _^_ 

l l (P L ,.,. t-h ( . (CIS) 



[p° - %-] + r( P ) 2 /4 V 2 

where r(p) = — 2Q : m^Sy i? (p) + Ef(p)), tanh = 1 — 2f(2p°), and ^ = eji with e^. 
|p| 2 /(2m*) and m* the non-relativistic effective mass. 
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